Spectral Curves and Whitham Equations in 
Isomonodromic Problems of Schlesinger Type 

Kanehisa Takasaki * 
Department of Fundamental Sciences, Kyoto University 
Yoshida, Sakyo-ku, Kyoto 606, Japan 
E-mail: takasaki@yukawa . kyoto-u .ac.jp 

Abstract 

It has been known since the beginning of this century that isomonodromic prob- 
lems — typically the Painleve transcendents — in a suitable asymptotic region 
look like a kind of "modulation" of isospectral problem. This connection between 
isomonodromic and isospectral problems is reconsidered here in the light of recent 
studies related to the Seiberg-Witten solutions of N = 2 supersymmetric gauge 
theories. A general machinary is illustrated in a typical isomonodromic problem, 
namely the Schlesinger equation, which is reformulated to include a small param- 
eter e. In the small-e limit, solutions of this isomonodromic problem are expected 
to behave as a slowly modulated finite-gap solution of an isospectral problem. The 
modulation is caused by slow deformations of the spectral curve of the finite-gap 
solution. A modulation equation of this slow dynamics is derived by a heuristic 
method. An inverse period map of Seiberg-Witten type turns out to give general 
solutions of this modulation equation. This construction of general solution also 
reveals the existence of deformations of Seiberg-Witten type on the same moduli 
space of spectral curves. A prepotential is also constructed in the same way as the 
prepotential of the Seiberg-Witten theory. 
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1 Introduction 



The notion of "isomonodromic deformations" was first discovered by R. Fuchs |T| in 1907 
as a new interpretation of the 6th Painleve equation (Pvi), and developed in diverse 
directions in the next decade. R. Fuchs considered a second order scalar ODE 

d 2 y 



d\ 2 



+ p(X)y = 



of Fuchsian type with four regular singularities, one of which is an apparent singularity. 
Gamier 0] generalized the work of R. Fuchs in two different forms. One generalization is to 
consider more than four regular singularities. This leasd to a multi-variable generalization 
of Pyi- The other is to include irregular singularities. The other five Painleve equations 
(Pi - Py) can be derived from this generalization. Schlesinger [H] obtained the so called 
Schlesinger equation from isomonodromic deformations of first order matrix ODE 

§ = M ^ 

where M(A) is an r x r matrix of the form 

N 4 . 

M(A)=E 



i=i 

Gamier noticed later 0] that the 2x2 Schlesinger equation is equivalent to his isomon- 
odromic deformations of 2nd order scalar Fuchsian equation. 

It is also in this decade that a link with "isospectral deformations" was uncovered. 
This is again due to Gamier || . He proposed an autonomous analogue of the Schlesinger 
equations, and pointed out that it can be integrated by Abelian functions. Remarkably, 
Gamier substantially arrived at the notion of isospectral deformations therein. Let us 
briefly review Garnier's discovery. Schlesinger's equation can be written 



9tj ti tj k(¥4) ^ 

This is a non-autonomous system, because the right hand side contains the independent 
variables t = {ti} explicitly. Garnier's proposal was to replace these tj's by constants, 
^ — > Q. The outcome is an autonomous system of the form 



Olj Ci Cj ht-f-;\ Cfc 
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As Gamier noticed, this gives isospectral deformations of M(A), namely, the characteristic 
polynomial det^M(A) — fil^j is invariant under the t flows. The algebraic curve defined 
by the characteristic equation 

det (M(A) -/i/) = 

in the (A, /i) plane (and its appropriate compactification) is nowadays called the "spectral 
curve." What Gamier did is to solve the above autonomous system in terms of Abelian 
functions on this algebraic curve. 

Isospectral problems of the same type were studied by Moser in the 70's. Moser 
proposed a matrix nonlinear system as a unified framework for a number of classically well 
known completely integrable dynamicals systems, such as the Neumann and Rosochatius 
systems, geodesic flows on an ellipsoid, etc. Moser's idea and Munford's related work 
were reformulated by the Montreal group (Adams, Harnad, Hurtubise and Previato) 
0, |9| and Beauville |H| to the isospectral problem of a rational matrix of the form M(A) 
(or U + M(A), where U is a constant matrix). Harnad and his collaborators later applied 
their method to isomonodromic problems |TT|, [Tj| [13], TA . 



Garnier's proposal, however, originally aimed at a quite different issue. He considered 
the autonomous system as a tool for studying asymptotic behavior of solutions of the 
Schlesinger equation in a neighborhood of singularities. A similar problem concerning the 



Painleve equations had been pursued by Boutroux Boutroux obtained an asymptotic 
expression of Painleve transcendents as a "modulated" elliptic function. Here "modula- 
tion" means that parameters of the elliptic function also depend (but "slowly") on the 
independent variable. In Garnier's program, Boutroux's elliptic curve is replaced by a 



more general algebraic curve. Flaschka and Newell |Tb| revisited Garnier's program in 
their study of isomonodromic and isospectral problems, and noted an important remark: 
They pointed out that a JWKB approximation converts the monodromy problem into a 
spectral problem. 

We reconsider this issue in the light of researches in the last ten years. Since the 
end of the 80 's, low dimensional string and topological field theories have provided new 
subjects of isomonodromic problems. A central subject is the so called "string quations" of 
two-dimensional quantum gravity JT7|. They are the first (or second) Painleve equation 
and its higher dimensional generalizations. These string equations have been studied 
from several different points of view. Among them, we are particularly interested in the 
approach [|T^, |T^, |2T|] from the "Whitham averaging method" (also referred to 



3 



as the "nonlinear JWKB method" |2I|). This method may be viewed as a modernized 



version of Boutroux's analysis. We now consider the Schlesinger equation in the same 
philosophy. 

The Whitham averaging method, or the nonlinear JWKB method, covers a wide area 
of modulational phenomena in nonlinear waves. The most relevant for our problem is the 
case of modulation of Abelian function solutions (which are usually called "finite-gap" or 
"quasi-periodic" solutions |25|) to a soliton equation. An "unmodulated" finite-band 
solution generally takes the form 



u o = Uo{£2Uiti | {In})- 



where Uq is an Abelian function, U^s are g-dimensional constant vectors (g is the genus of 
the spectral curve), and 7 n 's are other parameters of the solution. The parameters C/j's and 
J n 's are eventually determined by the spectral curve (as period integrals of meromorphic 
differentials). If the problem in question contains a small parameter e, one may consider 
a solution with the following asymptotic form 

u ~ u (l>i(?% I {/„}) 

i 

as e — > oo. The parameter C/j and I n now depend on the "slow variables" 

T = {Ti}, Ti = eU. 

This is the "modulation" of a finite-gap solution. The Whitham averaging method is a 
method to determine this slow dynamics in Tj in the form of differential equations. This 
kind of differential equations are generally called "modulation equations" (or "Whitham 
equations" ). In the case of finite-gap solutions, the modulation equation can be formulated 
as a dynamical system on the moduli space of spectral curves. 

These modulation equations of finite-gap solutions are known to possess a number of 



remarkable properties p6[ . It is Flaschka, Forest and McLaughlin [[27] who first pointed 
out that this type of modulation equations have a universal structure. They demonstrated, 
in the case of the KdV equation, that the modulation equation boils down to the universal 
form 

dTi 3 dTj ' 

where dfl^s are meromorphic differentials on the spectral curve. Krichever [^HJ and 
Dubrovin presented an abstract reformulation of this type of equations ("Whitham 
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hierarchies"), and constructed many special solutions (with applications to geometry and 
physics) . 

Our concern lies in special solutions of the universal Whitham hierarchy that represent 
slow dynamics of a spectral curve in isomonodromic problems. Remarkably, it seems likely 
that this class of solutions of the universal Whitham hierarchy are always characterized 
by a differential equation of the form 

d 

— — dS = dfli, dS = udX. 



This is indeed the case for the string equations fll8| , [19| , [20| , |2j| . In this paper, we shall 
derive a modulation equation of this form from the Schlesinger equation. 

We use a very heuristic method to derive the modulation equation. This heuristic 



method was developed in an attempt ]30| at an isomonodromic interpretation of integrable 
structures in supersymmetric gauge theories [31]]. The modulation equation turns out to 
possess almost the same properties as the so called "Seiberg-Witten solutions" of N = 2 
supersymmetric gauge theories. In particular, we introduce a period map of Seiberg- 
Witten map, and prove that the inverse period map solves the modulation equation. 
This also reveals the existence of another set of commuting flows on the moduli space of 
spectral curves. We also show that the notion of prepotential can be generalized to this 
case. These results will be strong evidence for the validity of the heuristic derivation. 

This paper is organized as follows. Sections 2 and 3 are of preliminary nature. In 
Section 2, we review basic properties of the Schlesinger equation. In Section 3, we consider 
the geometric structure of spectral curves along the lines of approach by the Montreal 
group and Beauville. Section 4, 5, and 6 are focussed on the derivation of the modulation 
equation. We begin with a reformulation of the Schlesinger equation in Section 4. The 
reformulated Schlesinger equation has a small parameter e. Garnier's autonomous system 
emerges in the limit of e — > 0. Our modulation equation is derived in Section 5, along 
with comments on other possible approaches. The structure of meromorphic differentials 
dQi, which are also basic constituents of our modulation equation, is specified in Section 
6. Section 7 is devoted to solving the modulation equation by the inverse period map. 
Section 8 deals with the notion of prepotential. We conclude this paper in Section 9. 
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2 Schlesinger Equation 



In this section, we review basic properties of the Schlesinger equation. For details and 
related topics, we refer to a series of papers by the Kyoto school (Jimbo, Miwa, Mori, 



Sato and Ueno) |3^, |33], |34 



2.1 Coadjoint Orbit and Hamiltonian Structure 

Let gl(r, C)^ denote a direct sum of N copies of gl(r, C). This is the space of iV-tuples 
(Ai, ■ ■ • , An) of r x r matrices. GL(r, C) acts on this space by the diagonal coadjoint 
action: Ai i— > gAig^ 1 . The Schlesinger equation can be written 



dtj 



(2.1) 



A i ,(l-5 ij )-^--5 ij j:-^ 

each coadjoint orbit Oi is left invariant under the t-flows. Thus the Schlesinger equation 
is actually a collection of non-autonomous dynamical systems on a direct product 0\ x 
■ ■ ■ x On of coadjoint orbits in gl(r, C). Usually, only semi-simple orbits are considered; 
such an orbit is labeled by the eigenvalues 6i a (a = l,---,r) of Ai. In other words, 
these eigenvalues (and, in general, the Jordan canonical form of A^s) are invariants of the 
Schlesinger equation. 

Actually, there are some extra invariants. They are the matrix elements of 

N 

^oo = — JZ^i; (2-2) 

i=l 

which are invariant under the Schlesinger equation: 

dA n 



du °- (2 - 3) 

This matrix, too, is usually assumed to be semi-simple, and it is customary to diagonalize 
this matrix in advance by a constant "gauge transformation" Aj t— > CAiC^ 1 . Thus only 
the eigenvalues (a = 1, • • • , r) of A^ are nontrivial invariants. 

Geometrically, this gauge-fixing may be interpreted as the Marsden-Weinstein con- 
struction of a "reduced phase space." The "unreduced phase space" is a coadjoint orbit 
Oi x ■ • • O n x Coo in the vector space gl(r, C) N+1 of (N + l)-tuples (Ai, ■ ■ ■ , A N , A^). In 
order to reproduce the Schlesinger equation, one has to impose the linear constraint 

N 

J2A + A oo = 0, (2.4) 
i=i 
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and "gauge away" redundant degrees of freedom by the action of GL(r, C). The left hand 
side of the linear constraint is essentially a moment map of this diagonal GL(r, C) action. 

Analytically, as we see show below, the coadjoint orbit invariants 6i a (i — 1, • • • , N, oo, 
a — 1, • • • , r) give local monodromy exponents of Schlesinger's monodromy problem. 

The coadjoint orbit structure leads to a Hamiltonian formalism of the Schlesinger 
equation. Let us introduce a Poisson structure on the vector space gl(r, C)^ by defining 
the Poisson bracket of matrix elements of A { = (A i>a p) a s: 

{Ai^ap, Aj pp} = 5ij(^—5fl p Ai yCta + 5 aa Ai p pj . (2-5) 

In each component of the direct sum, this is just the ordinary Kostant-Kirillov Poisson 
bracket. The Schlesinger equation can be written in the Hamiltonian form 

1 ={A j ,H i }, (2.6) 



dti 

where the Hamiltonians are given by 



and involutive, 



Hl = Res \ Tr M(A) 2 = £ Tr (^f\ , ( 2 .7) 



{Hi,Hj} = 0. (2.8) 



2.2 Isomonodromic Deformations 

The Schlesinger equation gives isomonodromic deformations of the first order ODE 

^ = M(X)Y (2.9) 

with the rational coefficient matrix 

N A. 

M(A)=E^ (2-10) 

i=l 1 

Note that A^ is the residue of M(A) at A = oo. 

Usually, this type of isomonodromic problems are considered under the following 

Assumption 

• The residue matrices Ai {% — 1, • • • , N, oo) are diagonalizable. 
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• The eigenvalues 9i a (a — 1, • • • , r) of each residue matrix Ai have no integer differ- 
ence, i.e., 9 ia — 9ip £ Z if a ^ f3. 

We assume them throughout this paper. These assumptions ensure that local solutions 
at the singular points \ = ti, ■ ■ ■ ,t]y, oo develop no logarithmic term (see below). 
The isomonodromic deformations are generated by the deformation equations 

dY Ai 
dU 



X-U 



-Y. 



(2.11) 



These deformation equations and the above first order ODE comprise an "auxiliary linear 
problem" of the Schlesinger equation. Its Frobenius integrability conditions can be written 
in the "zero-curvature form" 



9 



dtj X-U 



0, 



d_ Ai d_ 
dh X — ti dtj 



+ 



X - tj 



o, 



(2.12) 



and one can easily check that these zero-curvature equations are equivalent to the Schlesinger 
equation. 



2.3 Local Solutions at Singular Points and Tau Function 

Since A = t±, • • • ,tN and A = oo are regular singularities of the above first order ODE, 
one can construct a local solution of the following form at each of these singular points: 

• Local solution at X — tf. 

oo 

Yi — Yi-(X — Uf>, Yi = J2 y^n{X - U) n . (2.13) 

n=0 

• Local solution at A = oo: 

oo 

^oo = Yqq - A Yoq = ^2 ^oonA (2-14) 

n=0 

Here Y in are r x r matrices, the leading coefficients Y i0 and Y^q are invertible, and Oj 
and Oqo are diagonal matrices of local monodromy exponents. Inserting these expressions 
into the first order ODE gives the relations 

A = Y t0 e t Y^, i = l,---,7V,oo. (2.15) 

In particular, local monodromy exponents coincide with the eigenvalues of Af. 

e i = ti&g(e il ,--',e ir ). (2.i6) 
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It is not hard to check that these expressions of Ai and 0; are also consistent with the 
other equations of the auxiliary linear problem. 

The r function of the Schlesinger equation is defined in two equivalent ways. One way 
is to define logr as a potential of the Hamiltonians Hf 

N 

d\ogr = Y J H l dU- (2.17) 

i=i 

Another equivalent definition, which is more suited for generalization, is based on the 
equations 

^ = (2.18) 
The equivalence can be verified as follows: 

A { Aj 



Hi = V Tr 

A- 

= Res Tr - — l —M(X) 
A— tj A — tj 

= ResTr ^g^ ^y- 1 

= Tre.y^-'ya. 

The closedness of the 1-form J2HidU, or equivalently the integrability condition 

d Jk- d Jk (2 19) 

is ensured by the Schlesinger equation itself. 



3 Spectral Curve 

By "spectral curve," we mean the plane algebraic curve defined on the (A, //) plane by 

det(M(A) -///) =0 (3.1) 

and its suitable compactification. We first discuss its roles in isomonodromic and isospec- 
tral problems, then consider its geometric properties. 
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3.1 Spectral Curve in Isomonodromic Problem 



Isomonodromic deformations such as the Schlesinger equations are non-isospectral, namely, 
the spectral curve varies in deformations. 

Let us present an interesting formula (essentially due to Vereschagin |3l|) which show 
qualitatively that the characteristic polynomial of M(A) varies under isomonodromic de- 
formations. First, by the the well known identities of linear algebra, we have 

— logdet(M(A)-/i/) = — Trlog(M(A) - ///) 



By the first equation of the zero-curvature representation, dM(X)/dti can be rewritten 



dM{\) 
dU 



A: 



(A - 



A; 



x-u 



M(A) - ill 



The second term on the right hand side has no contribution, because of the identity 



Tr 



A; 



M(A) - fxl 



(M(A) 



0. 



.A - 1 

Thus, eventually, we obtain the formula 

9 log det (M(A) - fil) = Tr ^ (M(A) - /if) ~\ 



(3.2) 



This clearly shows that the t-dependence of the characteristic polynomial of M(A) is 
driven by the "anomalous" term Aj/(A — U) 2 . 



3.2 Spectral Curve in Isospectral Problem 

We now turn to Garnier's autonomous analogue of the Schlesinger equation. This equation 
has the following zero-curvature representation: 



8 A, 
\ — . 

dti X — c, 



M(A) 



0. 



A; 



d_ 

dti X — d ' dtj ' A — c,- 



d A, 
+ 3 



0. 



Repeating the same calculations as above, we now find that 

J-det(M(A)-/i/) =0, 



(3.3) 



(3.4) 



because there is no "anomalous" term like Ai/(X — ti) 2 . Thus one can confirm that Gar- 
nier's autonomous system is indeed an isospectral problem. An auxiliary linear problem 
is given by 



^ = M(A)V, |f 
dU 



A; 



A — c 



4- 



(3.5) 
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Here ip is understood to be a column vector. The first equation means that ip is an 
eigenvector of M(A) with eigenvalue /i. The second set of equations generate isospectral 
deformations. 

It is nowadays well known that this type of isospectral problems can be mapped to 
linear flows on the Jacobian variety of the spectral curve |24|, |25[ . Reproducing the original 
nonlinear problem is identical to Jacobi's inversion problem, and indeed solvable by theta 
functions (or Baker- Akhiezer functions). The solution thus obtained is written in terms 
of period integrals and Abelian functions. 

Although isomonodromic problems (in a generic case) cannot be solved in that way, 
the notion of spectral curve still plays a role in the study of Hamiltonian structures. 
Indeed, Harnad and Wisse [II, |I3J presented a construction of special Darboux coordinates 
( "spectral Darboux coordinates" ) using the language of spectral curves. The case of r = 2 
is particularly interesting, because this is the case where the 6th Painleve equation (Pvi) 
and Garnier's multi-variable version of Pyi (called the "Gamier system" by Okamoto 
||36|| ) emerge. The spectral Darboux coordinates in this case coincide with Okamoto's 
Darboux coordinates for the "Gamier system," which Okamoto discovered without using 
the notion of spectral curve. 



3.3 Geometry of Spectral Curves 

The structure of spectral curves of the above type has been elucidated in detail by the 
Montreal group M, |S| and Beauville [jTIJ. We present basic part of their results, which 
will be used in the subsequent sections. In the following, the poles of M(A) are written 
Cj rather than Tj. 

3.3.1 Spectral Curve on Plane. 

Let Co be the spectral curve on the (A, fi) plane: 

F(X, At) = det(M(A) - /il) = 0. (3.6) 

This becomes a ramified covering of the punctured Riemann sphere, 

7r : C -»• CP 1 \ {ci, c N , oo}, vr(A,/i) = A. (3.7) 

For a generic value of A, the inverse image 7r -1 (A) consists of r points (A, /j, a ) (a = 1, ■ ■ • , r). 
The /^-coordinates \i a of these points are eigenvalues of M(A). In a neighborhood of A = c«, 
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/i a 's behave as 

ji a = — — h non-singular, (3.8) 

A - a 

where 9i a are the eigenvalues of A^. Similarly, in a neighborhood of A = oo, 

» a = -e ooa \- 1 + 0(\- 2 ), (3.9) 

where 9ooa are the eigenvalues of = — Y.f=i M- (Of course, the numbering of /i a 's is 
meaningful only locally.) 

3.3.2 Compactification of Spectral Curve 

We now compactify C by adding several points over the punctures of the Riemann sphere. 
In a neighborhood of A = q, let us consider the following \x in place of //: 

N 

ji = f(\)li, /(A) = n(A-Q). (3.10) 

i=i 

In terms of the new coordinates (A, p), the equation of the spectral curve becomes 

F(X,jl) = det(/(A)M(A) - jll) = 0. (3.11) 

The inverse image 7r _1 (A) consists of r points (A, jj, a ) (r = 1, • • • , r) such that 

p. a = f'(.c i )e ia + 0(\-c i ) (3.12) 

as A — > q. Since the eigenvalues 9 ia (a = l,---,r) are pairwise distinct (recall the 
assumptions in Section 2), we add to C extra r points (A, ft,) = (q, f'(ci)9 ia ) to fill the 
holes above A = q. 

Similarly, in a neighborhood of A = oo, we use 

fx = A/i (3.13) 

in place of /i. 7r _1 (A) now consists of the r points (A, j2 a ) (r = 1, . . . , r) such that 

ft a = -9 ooa + 0(\- 1 ) (3.14) 

as A — > oo. The eigenvalues 9^ (r = 1, • • • ,r), too, are pairwise distinct. Therefore we 
add to Co the r points (A,/i) = (oo, — #ooa) to fill the holes above A = oo. 

Thus, by adding altogether rN + r points to C , we obtain a compactification C of C . 
The covering map 7r uniquely extends to C, and gives a ramified covering 7r : C — > CP 1 
of the Riemann sphere. 
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3.3.3 Genus of Compactified Spectral Curve 

If there are multiple eigenvalues of the compactified spectral curve C has singular 
points (nodes) over A = q. In that case, one has to take a desingularization p : C — > C 
for further consideration. Fortunately, this does not occur in our case because of the 
assumptions introduced in Section 2. One can show, by a standard method, that C has 
the genus 

g = -(r- l)(rN-r-2). (3.15) 

3.3.4 Structure of Characteristic Polynomial 

In order to examine the structure of the characteristic polynomial of M(A) (or, rather, 
/(A)M(A)), we take the rational matrix 

N AO 

M°W = ^, A° = diag(^,---,^), (3-16) 
l=l A - q 

as a reference point on the coadjoint orbit that M(A) belongs to. Now compare the 
characteristic polynomials F(\,p) and F°(X,p) of /(A)M(A) and /(A)M°(A). Note, first, 
that F(\,fi) — F°(A,/i) vanishes at A = q. This will be obvious from the following 
expression of these polynomials as A — > q: 

F(A, //) = det(/'(q)A - fil) + 0(A - q), 
F°(A,/i) = det(/'(q)A° - jil) + 0(A - q). 

(The two determinants on the right hand side are equal, because Ax and A\ are on the 
same coadjoint orbit.) Therefore F(X,p) — F°(A,/i) is divisible by /(A). Furthermore, 
these characteristic polynomials have the following expansion in powers of //: 

F(A, /2) = (-^) r + Tr M(A)(-/i) r ~ 1 + • • • , 
F°(A, /i) = (-/^ + Tr M°(A)(-/ir 1 + • • • . 

Since the first two terms on the right hand side are equal, respectively, the difference of the 
two polynomials contains no terms proportional to fi r and /i r_1 . From these observations, 
one can conclude that the difference of the characteristic functions can be written 

F(A,/i)-F(A,/i) = /(A)^MA)/i r - £ , 
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where p^(A)'s are polynomial functions of A. A simple power counting argument (assigning 
weight 1 to A and weight TV — 1 to p) shows that the degree of pe(X) does not exceed the 
positive integer 

8 e = (N-l)e-N. (3.17) 

Therefore pe(X) can be written 

h 

m=0 

The leading coefficient of this polynomial is a function of ^ a 's only, because it can be 
writen 

h sa = (-lY^eiA^) - a e (A°J), (3.18) 

where a^A^) and ae(A^ X3 ) are the £-th elementary symmetric function of A^ and A 1 ^ = 
— Y7i=i A®. Thus, eventually, we arrive at the following expression of F(X,jl): 

F(A,/2) = F°(X,jl) + £ £ h m iX m jjL r ~ e . (3.19) 

1=2 m=0 

3.3.5 Parameters of Spectral Curve 

Apart from the leading coefficients, the coefficients h m £ (£ — 2, • • • , r, m — 0, • • • , Se — 1) 
of p^(A)'s give arbitrary parameters ("moduli") of the spectral curve. Their total number 
coincides with the genus of the spectral curve: 

±5 i = j:((N-l)£-N)=g. (3.20) 

1=2 1=2 

In the isospectral problem, these parameters h m £ are constants of motion (Hamiltonians 
of commuting isospectral flows). They should not be confused with the Hamiltonians Hi 
of the isomonodromic deformations. 

Actually, the characteristic polynomial has yet another set of parameters — the posi- 
tion Q of poles of M(X). They play the role of deformation variables in isomonodromic 
problem. 

In summary, the spectral curve has three distinct sets of parameters: 

• Position of poles q (i — 1, • • • , N). 

• Coadjoint orbit invariants 9i a (i — 1, • • • , N, oo). 

• Isospectral invariants h m i (£ — 2, • • • , r, m — 0, . . . , 5g — 1). 
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4 Schlesinger Equation with Small Parameter 



4.1 Reformulation Including Small Parameter 



The first step towards the derivation of our modulation equation is to reformulate the 
Schlesinger equation by the following substitution rule: 



d_ 



a 



d_ 

dX 



d 



A; 



A; 



(4.1) 



&Ti OX «9A' X-U X-Tt 
Note, in particular, that the deformation variables are renamed as ti — > Tj. Tj's will play 
the role of "slow variables." The reformulated Schlesinger equation reads: 



:i - St 



E 



[Ai, Ak 



rji rji / j rp rp 

-Li - J-j k ^ {) J-i - J-k 



An auxiliary linear problem is given by 



dY dY 



A; 



-Y. 



(4.2) 



(4.3) 



A - Tj 

The above e-dependent Schlesinger equation can be indeed reproduced from the the Frobe- 
nius integrability conditions 



_d_ 



d 

e m + x 



d 



A, 



0. 



0. 



(4.4) 



Ti dTj X - Tj 
A few comments on the above reformulation will be in order: 

(i) Reformulating the Schlesinger equation as above is inspired by the work of Verescha- 
gin |3"5| , who considered all the six Painleve equations (Pi-Pyi) in such an e-dependent 
form. This is different from the way Boutroux JTj| and Gamier [5[] derived an Abelian 
function approximation, but, as Vereschagin stresses, they are asymptotically related. 

(ii) The string equations of two-dimensional quantum gravity, which are Pi, Pn and 
their higher order generalization, contain such a small parameter from the beginning. The 
small parameter i interpreted as "string coupling constant" IT 7 



4.2 Isospectral Deformations as First Approximation 

One can now see, at least intuitively, that an isospectral problem emerges from the above 
e-dependent Schlesinger equation. Let us introduce the "fast variables" 

U = e- 1 ^, (4.5) 
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and rewrite the above equation as 

dAi [Ai,Aj] ^ ( . „s 

— = (1 - 5ij) — - d y 2^ _ . (4.6) 

Suppose we now observe this system m the scale of the fast variables U. In this scale, 
Tj's may be treated as being approximately constant, because a finite displacement in tj's 
corresponds to a small (i.e., 0(e)) displacement in Tj's. If Tj's were true constants, the 
above equation would be exactly Garnier's autonomous system. Thus, observed in the 
scale of the fast variables, our e-dependent Schlesinger equation looks approximately like 
an isospectral problem. 

Of course, this is no more than an approximation. In fact, Tj's are not constant, but 
vary slowly in the order of 0(e). Accordingly, the spectral curve, too, deforms slowly 
because the defining equation of the spectral curve contains Tj's explicitly. A precise de- 
scription of the isospectral approximation has to take into account such slow deformations 
of the spectral curve. 

It is, however, not only Tj's that vary; the isospectral invariants h m £ also change 
values slowly. They are both responsible for deformations of the spectral curve. In order 
to see how this occurs, recall the derivative formula of the characteristic polynomial under 
isomonodromic deformations (Section 3.1). In the present setting, this formula takes the 
form 

— logdet(M(A) -pi) = Tr - - ' (iVf(A) - ///)". (4.7) 

The right hand side plays the role of a "driving force" for slow deformations of the 
characteristic polynomial. More precisely, in order to extract a true driving force of slow 
deformations, one has to take an average of the right hand side over the quasi-periodic 
motion in the fast variables tj. This is a central idea of Vereschagin's averaging method 



Meanwhile, the coadjoint orbit invariants 6i a remain constant, because they are also 
invariant of isomonodromic deformations. 

4.3 Concept of Multiscale Analysis 

Although frequently lacking mathematical rigor, "multiscale analysis" is widely accepted 
in applied mathematics as a powerful tool for dealing with this kind of problems |37] . 



The Whitham averaging method, too, has been developed in the framework of multiscale 
analysis. The nonlinear JWKB method of Dobrokhotov and Maslov is an attempt 
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at a rigorous reformulation of multiscale analysis. Let us show how our problem may be 
formulated in the language of multiscale analysis. 

A key of multiscale analysis is to treat the fast variables t = {U} and the slow variables 
T = {Tj} as independent variables. A^s are thus assumed to be a function of t and T 
(and the small parameter e), 

A i = A i (t,T,e). (4.8) 

The relation tj = e _1 Tj is imposed only in the differential equation in question. (For a 

more precise description, one may introduce a series of "slower" variables t^ 2 \ ■ ■ 
related with t fc 's as 

tf ] = e% (4.9) 

For our purpose, only the first two scales are sufficient.) 

In this multiscale ansatz, derivative terms in the equation are given by a sum of 
t-derivatives and T-derivatives. In terms of differential operators, this amounts to substi- 
tuting 

d d d , 

t wrm, +t m- (410) 

We now assume an asymptotic expansion of the form 

A t (t, T, e) = Af\t, T) + Af\t, T)e + ■ ■ ■ , (4.11) 

and plug all these stuff into the Schlesinger equation. From each order of e, we obtain a 
differential equation for the coefficients of the above expansion. 
The lowest order equation is give by 

a 4(0) [4(0) 4(0)1 r 4 (0) .(0)1 

^- = (1 - iJ-^A - E l 4fr^- (4.12) 

If we consider Tj = q, this is nothing but Gamier 's autonomous system. Note that 
this derivation of Garnier's isospectral problem is more understandable than the intuitive 
derivation in the last subsection. In the derivation of the last subsection, t and T were not 
independent and constrained by the relation tj = e _1 Tj, thereby we had to say that T is 
"approximately constant"; in the setting of multiscale analysis, t and T are independent. 
This shows a conceptual advantage of multiscale analysis. 

The lowest order equation, however, carries no information on the T-dependence of 
A^'s, which are to be determined by the next order equation. The next order equation 
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is given by 



(1 — 8n i 

dU 3 Ti-T, 



J 1 3 

^ _(o) , r t"»i 



+ (terms containing A^'s and their T-derivatives only). (4.13) 

This equation contains T-derivatives of Af^ as well as t-derivatives of A± . A standard 
prescription of multiscale analysis is to eliminate the latter by "averaging over the t 
space." This procedure usually takes the form of an "integrability condition" for the above 
equation to have a solution A^p under suitable a boundary condition (e.g., requiring the 
same quasi-periodicity as A^'s possess). 

One thus obtains a differential equation (in T) for t-averaged functionals of Af^s. 
This is one of various possible expressions ("modulation equations") of modulational 



dynamics. As Whitham pointed out such a modulation equation frequently appears 



in the form of "averaged conservation laws." Having this fact in mind, Flaschka, Forest 



and McLaughlin |27j considered averaged conservation lows of the KdV equation, and 



derived their compact expression of this problem. 

4.4 Our approach to Modulation Equation 

There are many technical difficulties in deriving a modulation equation from the Schlesinger 
equation along the line presented above. A main obstacle is the fact that the spectral curve 
is no longer hyperelliptic for r > 2. Most attempts in the literature, including the work of 
Flaschka, Forest and McLaughlin, have been limited to hyperelliptic spectral curves. This 
considerably reduces the obstacles (though a complete treatment of the problem is still 
by no means an easy task). The spectral curve of the Schlesinger equation is hyperelliptic 
only if r = 2. In the general case, the averaging method is inevitably confronted with 
delicate problems of geometry of spectral curves. (In this respect, Krichever's averaging 
method [|38| for general spectral curves is quite remarkable.) 



Since our main concern is the structure of the modulation equation rather than the 
averaging method itself, we now bypass these delicate issues by a very heuristic argument. 

5 Modulation Equation 
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5.1 Multiscale Analysis of Auxiliary Linear Problem 



The first step of our heuristic argument is to apply the concept of multiscale analysis to 
the auxiliary linear problem of the Schlesinger equation. Besides the multiscale expression 
of AiS, we now assume the following ansatz to Y (which is now understood to be vector- 
valued): 

Y = (0 (o) (t, T, A) + (1) (t, T, A)e + • • •) exp e^Sp, A). (5.1) 

Here (j>^(t, T, A) are vector- valued function and S(T, A) a scalar- valued function. This is a 
kind of JWKB ansatz (inspired by the work of Flaschka and Newell |TjJ and Novikov [|TB|). 
Note that the "phase function" S is independent of t. We can now write the auxiliary 
linear problem in the following multiscale form: 

dY 



dX 



' d_ _d_ 



Y 



M(A)y, 

A: 



X — Ti 



-Y. 



(5.2) 
(5.3) 



The leading order equation should reproduce the auxiliary linear problem of the 
isospectral problem. Let us confirm that this is indeed the case. The leading order 
equation is given by 

dS{X) 



dX 



-0(°)(A) =MW(A)0(°)(A), 



(A) , ^(A) 0(O)(A) 



dti 



+ 



A 



(o) 



X - Ti 



where 



We now define 



N AO) 
i=l A - 1 i 



dS(X) 

~dx~ 



^ = 0(°)(A) exp^t, 



N 



i=l 



dS{X) 



(5.4) 
(5.5) 

(5.6) 
(5.7) 



In terms of these quantities, the leading order equation of multiscale expansion can be 
rewritten 



^ = M(°)(A)^, ^ 



A- 



(0) 



-1>. 



(5i 



X — Ti 

This is exactly the auxiliary linear problem of the isospectral problem that we derived in 
Section 4! 
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5.2 Matching with Baker- Akhiezer Function 

The next step is the most crucial part of our heuristics. We now compare the above ip 
with the Baker- Akhiezer function in the ordinary algebro-geometric approach |24|, 25 1 to 
finite-band solutions of soliton equations. 

It is well known that the solution of the auxiliary linear problem for a finite-band 
solution can be constructed as a (scalar- or vector-valued) Baker-Akhiezer function. In 
the present setting, such a Baker-Akhiezer function can be written 

N 

ip = 4>exp^2tiQi. (5.9) 

Here is a vector-valued function on the spectral curve, also depending on t; each entry 
is a combination of theta functions. We do not specify its detailed structure, because it 
is irrelevant in the following consideration. Meanwhile, Qi is the primitive function of a 
meromorphic differential dQi on the spectral curve, 

Qi = / dfli. (5.10) 



Note that these quantities also depend on T through the T-dependence of the spectral 
curve. Another important remark is that such an expression of the Baker-Akhiezer func- 
tion is available only after selecting a "symplectic homology basis" of the spectral curve, 
i.e., cycles Ai, Bj (I = 1, • • • , g) with intersection numbers Aj ■ Aj = Bj ■ Bj = and 
Aj ■ Bj = 5ij. (This issue will be discussed in detail in Section 6.) 

We now assume that, for a suitable symplectic homology basis, this Baker-Akhiezer 
function coincides with the ip derived from the leading order equation of multiscale expan- 
sion. More precisely, we require their "amplitude part" and "phase part," respectively, to 
coincide: 

(o) = 0, || = n*. (5.ii) 

Actually, the first relation may be relaxed as 

<f>W = <f>h(T,\), (5.12) 

where h(T, A) is a scalar function independent of t; this takes into account the obvious 
symmetry 

0(0 ^(f)( )h(T,X) (5.13) 
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of the lowest order equation of multiscale expansion. Such a factor h(T, A) might be 
necessary to proceed to the next order approximation of multiscale analysis. We shall not 
go into this issue here. 

Remark. We take this opportunity to correct an error in the previous publication 
3~0fl . Substantially the same JWKB ansatz also assumed therein, but with an extra 



factor of the form (d 2 S/dX 2 ) 1 ^ 2 — see Eqs. (13) and (21) therein. This factor itself is 
nonsense. Fortunately, this does not affect the leading order of multiscale analysis. As 
mentioned above, however, something like this factor will be necessary in the next order. 
The arbitrary function h stands for such a compensating factor. 

5.3 Modulation Equation 

Thus, under several assumptions, we have been able to derive a series of relations that 
link the isomonodromic and isospectral problems. In particular, the following equations 
are obtained: 

dx=^ m=^ (5 - 14) 

We propose these equations as the modulation equation that governs slow dynamics of 
the spectral curve 

det (M (0) (A) -fil) = 0. (5.15) 

Remember that one has to select a symplectic homology basis Ai,Bi (I = 1, • • • ,g) 
in order to derive these equations. Actually, selecting a proper homology basis is a non- 
trivial problem. We discuss this issue in more detail in Section 6, along with a precise 
characterization of the meromorphic differentials dQ{. 

5.4 Relation to Generic Whitham Equation 

As a piece of evidence that our modulation equation is a reasonable one, we now show 
that solutions of this equation give a special subfamily of generic Whitham deformations. 
The first equation of the modulation equation can be rewritten 

dS = —dX = fidX. (5.16) 

CM 

In other words, S, like Q^s, is the primitive function of a meromorphic differential dS on 
the spectral curve: 

S= / fjdX. (5.17) 



21 



The second part of the modulation equation implies the equations 



-^-dS = dVL i . (5.1* 



Accordingly, the generic Whitham equation 



dflj = -=-dSl 5 (5.19) 



dTi J dT 3 



follows immediately. (Here, as usual, A is understood to be constant under the T- 
derivation. Geometrically, this should be treated as a connection. See the paper of 
Krichever and Phong [j39|j .) 



Thus, our modulation equation turns out to yield a special subfamily of generic 
Whitham deformations. We shall see in Section 8 that an inverse period map yields 
generals solution of this equation. 

5.5 Possible Approach from Averaging Method 

Our heuristic argument has to be justified in a more rigorous form, or at least cross- 
checked by some other method. A possible check will be to derive the above equations by 
an averaging argument. 

Krichever's averaging method J38| (see also the papers by Fucito et al. |20| Carroll 



and Chang |4(|) seems to be particularly promising. This method employs the so called 
"dual Baker- Akhiezer function" if)*, i.e., a solution of the dual auxiliary linear problem 

l*l> m = PM®(\), -^-=^^ ¥ (5.20) 

along with ip. One can derive, for instance, the following formula evaluating T-dependence 
of the eigenvalue /x of M^>: 

a„ {4.'A«Xx-T t) -H) 

m~ < n- > ' ' 

Here < • ■ ■ > means the average over the fast variables t. This formula can be readily 
translated to the language of dS (though we have been unable to calculate the average 
and to identify the resulting equation with the above equation for dS). Note that the 
term A^(\ — Tj)~ 2 in the above formula is the same as the "driving force" term that we 
encountered in the derivative formula of the characteristic polynomial of M(A). 
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6 Structure of Meromorphic Differentials 

In this and subsequent sections, we omit the superfix and write Af\ M^(A), etc. 
as Ai, M(A), • ■ •. This is just for simplifying notations. 

6.1 General Remarks 

The vector-valued Baker- Akhiezer function if) is single valued on the spectral curve, mero- 
morphic outside the points of 7r _1 ({Ti, • • • , T/v}), and has essential singularities of a par- 
ticular exponential form at these exceptional points. The entries of <j) have a common pole 
divisor D, part of which may overlap with the essential singularities. (In the following, we 
consider a generic case where this overlapping does not occur.) The essential singularities 
at points in 7r _1 (Tj) are generated by poles of the meromorphic differential dQi. 

The meromorphic differentials dfli have to be selected for the essential singularities 
of ip to match the auxiliary linear problem. This yields conditions on the poles of the 
meromorphic differentials. Our task in the following is to specify those conditions. 

In order to consider this problem, it is convenient to build a matrix solution \1> of the 
auxiliary linear problem from the vector- valued Baker- Akhiezer function ip. Let P a (\) 
(a = l,---,r), denote the points of 7r _1 (A). (Of course, also here, such numbering is 
meaningful only locally.) ^ is given by 

*= (^(Pi(A)),...,^(P r (A))). (6.1) 

This matrix has the following factorized form: 

^ = $expdiag(^t i a(Pi(A)),---,E^(Pr(A))). (6.2) 

$ = $(A) is an r x r matrix originating in the "amplitude part" 0, and invertible at a 
generic point. 

6.2 Conditions on Meromorphic Differentials 
6.2.1 Poles of dVLi 

All necessary information on poles of dVti can be derived from the following relation, which 
is obtained by inserting the above expression of ^ into the auxiliary linear problem: 

l^- 1 + $ diag(^(Pi(A)), • • • , SkiPrW))*- 1 = "jAv ( 6 - 3 ) 



23 



We now compare the singular part of both hand side at A = Tj. Since the first term 
on the left hand side is non-singular, it turns out that f2j has a pole of first order at each 
point of 7T _1 (Tj), and the residue is equal to (—1) times an eigenvalues of Since the 
eigenvalues of are 9 ia , this implies that (after suitably renumbering the eigenvalues) 
Ai can be written 

A = $(T,) diag(0 a , • • • , 9 ir )<f>(T i )-\ (6.4) 

and that 

Q 

fij = — %OL + non-singular (6.5) 
A — Tj 

in a neighborhood of P a (Tj). 

In particular, the meromorphic differential dfli behaves as 

9- 

dQi = — — ia 2 d\ + non-singular (6.6) 
(A - Ti) 

in the same neighborhood of P a (Tj). Therefore, if 9 ia ^ 0, P a (T;) is a pole of dVli of 
second order. These are all poles that dVti is required to have. 

6.2.2 Period Integrals of dVli 

The above conditions on poles determine the meromorphic differential dVti up to a dif- 
ference of holomorphic differential. We now select a symplectic homology basis Aj, Bj 
(I = 1, . . . , g), and put the standard normalization condition 

£ dft i = 0, I = l,...,g. (6.7) 

dfli is thus uniquely determined. 



6.3 Back to Modulation Equation 

This is the end of the precise description of our modulation equation. Let us finally 
reconfirm the roles of the three sets of parameters 9i a , Ti and h m g in the defining equation 
of the spectral curve: 

• The coadjoint orbit invariants 9i a are constant parameters. 

• The positions of poles Tj are "time variables." 

• The isospectral invariants h m £ are "dynamical variables." 
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6.4 Remarks on Symplectic Homology Basis 

The choice of symplectic homology basis is a non-trivial problem. The work of Flaschka, 
Forest and McLaughlin |27j provides a typical example for considering this problem. In 
their work, the meromorphic differentials dQi are normalized along the so called "/i- 
cycles." These cycles are homologous to trajectories of the "auxiliary spectrum," which 
are nothing but the degree g divisor Sf = i(Aj, Q) of zeros of a Baker- Akhiezer function on 
the KdV hyperelliptic spectral curve 

2g+l 

c 2 =n^-^)- (e.s) 

1=1 

The Whitham averaging over the quasi-periodic motion on the Jacobian variety eventually 
boils down to period integrals along these trajectories. As this typical example shows 
clearly, the normalization condition of dTti is by no means arbitrary, but determined by 
the geometric structure of motion of a divisor on the spectral curve. (In fact, as Ercolani, 
Forest and McLaughlin noted JIT], this issue is already considerably delicate in the case 
of the sine-Gordon equation.) 

Meanwhile, our formulation of the modulation equation itself works for any choice of 
the symplectic homology basis. Furthermore, the results of the subsequent sections also 
hold irrespective of the choice of the symplectic homology basis. Presumably, it is only 
a subset of solutions of our modulation equation that actually correspond to the true 
modulational description of isomonodromic problems. This issue, too, forms part of the 
hard analytical problems that we do not pursue in this paper. 



7 Solutions of Modulation Equation 

7.1 Period Integrals a la Seiberg-Witten 

In the following, let h = {hj | I — 1, • • • ,g} denote the g-tuple of isospectral invariants 
h m £ (£ — 2, • • • , r, m — 0, ■ ■ ■ , 8i — 1) (see Section 3). With the other parameters 9 ia and 
Tj being fixed, the spectral curve forms a ^-dimensional deformation family parametrized 
by these isospectral invariants. This is the same situation as the Seiberg-Witten solution 



and its various generalizations |31[| . One can indeed derive the following analogous results. 



A set of (local) coordinates on the ^-dimensional moduli space of spectral curves 
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are given by the period integrals 

ai = i dS, I = 1, • • • ,g. (7.1) 

J At 

Here (and in the following), Ai, Bj (J = 1, • • • , g) are the same symplectic homology 
basis as used in the definition of the modulation equation. 

Another set of (local) coordinates are given by the dual period integrals 

h=l dS, 1=1,- .-,0. (7.2) 
In the Seiberg-Witten theory, they are denoted by af. 

There is a (locally defined) function T = J 7 (a) ( "prepotential" ) of a = {a/} such 
that 

The second derivatives of T coincide with matrix elements of the period matrix, 

jr4- = r u=f d ^ ( 7 - 4 ) 

oajoaj jb 1 

where dui (I — 1, - • • , g) are a basis of holomorphic differentials uniquely determined 
by the normalization condition 

duj = 6 U . (7.5) 



A 



In a sense, the rest of this paper is devoted to verifying these results in an T-dependent 
form. 

The goal of this section is to show that the inverse period map a i— > h solves our 
deformation equation. We first establish the invertibility of the period map h i— > a. 
The inverse map a i— > h then turns out to satisfy a deformation equation of Seiberg- 
Witten type (with respect to o^'s), as well as our modulation equation (with respect to 
Tj's). Most results and proofs are a rather straight forward generalization of those of the 
Seiberg-Witten solution 

The notion of prepotential will be discussed in detail in Section 8. 
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7.2 Invertibility of Period Map 

This subsection is organized as follows. Firstly, we construct a basis dui (I = 1, • ■ ■ , g) of 
holomorphic differentials as derivatives of dS with respect to hi's. Secondly, we examine 
the linear relations between this basis and the normalized basis duj (I — 1, • • • , g). As a 
corollary, (local) invertibility of the period map h —>■ a follows. 

7.2.1 Basis of Holomorphic Differentials 

We here consider /x to be a (multivalued) function fi = /i(T, h, A) of T, h and A, and 
differentiate dS(T, h, A) = fi(T, h, X)d\ by hj. This gives the differential 

do, = » dS = MLhR dX . (7 . 6) 

ohi ohj 

We now show that this gives a basis of holomorphic differentials on the spectral curve. 
The following reasoning is borrowed from the work of Adams, Harnad and Hurtubise |§ . 

This definition of dujj can be cast into a more tractable form as follows. First rewrite 
it in terms of jl = /(A)/x (see Section 3): 

duJl = 77TY ai. dX - 

/(A) dh T 

Differentiating the equation 

F(\,jl) = det(/(A)M(A) - yl) = 

of the spectral curve in the (A, p) coordinates gives 

djl dF(X } jl)/dh I 
dh~i ~ ~ dF{X,p)/djl ' 

Now recall that hi is just an abbreviation of one of hmi's in the following formula of 
Section 3: ^ 

F(X, fi) = F°(A, fl) + /(A) £ J2 h m iX m fi r - £ . 

i=2 m=0 



Therefore, for hj = h m £, 

dF(X,fl) 



/(a)aV 



dh m £ 

Thus, eventually, we obtain the following expression of dui = duj m e: 



X^fL T ~ m 

duj m i = = dX. (7.7) 

dF{X,fi)/df l V ; 
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The last expression may be viewed as the Poincare residue of a meromorphic 2-form 
on the (A, ft) plane with pole divisor along the spectral curve. Since 2 < £ < r and 
< m < 5g — 1, a standard argument leads to the following result: 

Proposition 1 The differentials dui (I = 1, • • • ,g), form a basis of holomorphic differ- 
entials on the spectral curve. 

7.2.2 Linear Relations between Two Bases 

Since dwi and dui both give a basis of holomorphic differentials, they should be linked 
by an invertible linear transformation: 

9 

duj 1 = Aijduj. (7.8) 
j=i 

The matrix A = (Au) of the coefficients is invertible. 
In fact, the matrix A is given by the period integrals 



A 



u 



I dw I: (7.9) 
JAj 



as one can readily see by integrating the above linear relation of dui and duoj along Aj. 
Similarly, integrating along Bj yields the matrix relation 

B = AT, (7.10) 

where the matrix elements of B = (Bu) are given by 



B L j= i dwj. (7.11) 
JBj 

In particular, the period matrix T can be written 

T = A- 1 B. (7.12) 

7.2.3 Invertibility of Period Map 

We now show the following result, which shows that the matrix is in fact the Jacobian 
matrix of the period map h t— > a. 

Proposition 2 

(7,3) 
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Proof, a j is defined by the period integral 

a j = f dS. 



A 



J 



Its /^/-derivative can be calculates as follows: 

daj f d 



dS = <f> duji = Ai 



j 



dhi Jaj dhi Ja 
Q.E.D. 

Since A is an invertible matrix, we have: 
Corollary 1 The period map h i— > a is (locally) invertible. 

7.3 Solving Deformation Equations by Inverse Period Map 

We now take into account the variables T = {Tj}, and consider the i-th coordinate of the 
inverse map o h ft as a function of T and a: 

/i/ = /i/(T,a). (7.14) 

(In order to obtain these functions, one has to solve the defining equation of the a/s for 
the parameters hi. This is by no means an easy task.) Accordingly, dS = fxd\ becomes 
a meromorphic differential depending on the parameters (T,a): 

dS = dS{T,a). (7.15) 

We can now prove that dS satisfies the modulation equation with respect to T: 

Proposition 3 dS satisfies the modulation equation 

d 



-dS = da t . 



Proof . Let us tentatively define 



dVti = J^-dS = ^-d\. 
dTi dTi 

We show that this differential satisfies all conditions that characterize dQi. By the unique- 
ness, then, dVt-i coincides with dQi, and the proof is completed. 

(i) Locations of poles of dVti. Differentiating the equation of the spectral curve 



F(X, aO = det(M (A) - fil) = 
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gives the relation 

dfi dF(\,n)/dTj 

df~ ~ ~ dF(\,fj,)/dfi' 
Therefore, apart from the poles of fi, the derivative dfi/dTi can have poles at the zeros of 
dF/dfi (i.e., at ramified points of ir). Let e be the ramification index at a ramified point. 
dfj,/dTi has a pole of order at most e — 1 there. This pole, however, is canceled by zeros of 
dX because dX has a zero of order e — 1 at the same point. Meanwhile, in a neighborhood 
of each point of 7r _1 (oo), 



d/j, 

m 



(^ a A- 2 + 0(A- 2 ))rfA, 



so that the derivative d/i/dTi cancels the second order poles of dX at these points. Thus, 
poles of dQi are limited to points in 7r _1 ({T 1 , • • • , T/v}). 

(ii) Singular behavior of dQi at poles. We examine the singular behavior of dQi at 
each point of n~ l {Tj) = {Pji, • • • , Pj r }- Recall that, in a neighborhood of Pj a , /i behaves 
as 

\i = — - — ^7=- + non-singular. 

A - Tj 

Therefore 

dfi I non-singular (j ^ i) 

dTi | ~ ( A -'r f + non - sin g ular U = *) 

Consequently, dfli has poles only at the points over A = Tj, and exhibits there the same 
singular behavior as dQi does. 

(iii) Period Integrals of dQi. Since the deformations leave a/'s invariant, we have 

- ^FdS=±f dS = ^ = 0. 
Aj dTi dTi J Al dTi 

Thus dQi satisfies the normalization condition of periods, too. Q.E.D. 

Similarly, dS turns out to satisfy a deformation equation of Seiberg-Witten type with 
respect to a/s: 

Proposition 4 dS satisfies the deformation equation 

d 

dS = du) i, I — 1, • • • , g, 



daj 

of Seiberg- Witten type. 
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Proof. This is just a consequence of the chain rule of differentiation: 

——dS = V ^7— dS = y2{A~ l ) u dujj = du^ 
oai j oaj ohj j 

Q.E.D. 

Thus, we have been able to show that the inverse period map a ^ h solves the coupled 
deformation equations 

d d 
— —dS = dtti, ^r—dS = dujj. (7-16) 
oTi oaj 

(Reconfirm, once again, that the differentiation is understood to leave A constant, i.e., 
OX/dTi = dX/dd! = 0.) 



7.4 Remarks on Structure of dS 

It seems remarkable that although our meromorphic differential dS has simple poles with 
non-zero residues, their derivatives in deformation variables are holomorphic differentials 
or meromorphic differentials of second kind (i.e., meromorphic differentials with higher 
order poles and no residue). This is reminiscent of the structure of the Seiberg-Witten 



meromorphic differential in supersymmetric gauge theories coupled with matter fields |31 

Because of this fact, dS cannot be written as a linear combination of dwi and dQf, we 
need extra meromorphic differentials of the third kind, say dU, with the same singularity 
structure as dS but with vanishing ^/-cycles. dS can be written 

9 

dS = Y,aiduJi + dU. (7.17) 
1=1 

Of course dU can be further decomposed into a sum of more elementary meromorphic 
differentials, but can never be a linear combination of dQi's. 



8 Prepotential 

8.1 Definition of Prepotential 

The notion of prepotential has been formulated in a quite general framework by Krichever 
|28fl and Dubrovin [^3]. Following their formalism we can now define a prepotential 
T = J-(T, a) by the following differential equation. 
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Proposition 5 The following differential equations is integrable in the sense of Frobenius: 

^-=i dS, % = V Res -r^-dS. (8.1) 
d ai J Bl dTi £i Pic, A - Ti K J 

Here Pi a , a — 1, • • • ,r, denote the points of n~ l {Ti) such that 

q 

H = ta + non-singular (8.2) 
A — Ti 

in a neighborhood of P ia . 

8.2 Equivalent Definition 

Before proving the above result, we note here that the defining equations of T can be 
rewritten in the following more compact form: 

, rr 

The equivalence of the first equation is obvious from the definition of bj. Let us verify 
the equivalence of the second equation in some detail. We first rewrite the definition of 
Hi into a contour integral 

Hi = -== / - Tr M(Xfd\ 

along a sufficiently small circle (of radius S) around X — Ti. Let fi a , a = l,---,r, 

denote the eigenvalues of M(A); they correspond to the r sheets of the spectral curve in 

a neighborhood of 7r _1 (Tj). Expressing the trace of M(A) 2 in terms of these eigenvalues, 

we can rewrite the above integral formula as: 

1 r r u 2 

Hi 1 ^ c 



/|A-Ti|=tf^i 2 



27TV-1 J\\-Ti\=S 

= t^=f —dX. 
27r v /= T /c ia 2 

Now recall that /i has a Laurent expansion of the form 

q 

V = - + Ciafi + c ia ,i(A — Ti) H 

in a neighborhood of Pj Q . From this fact, we can easily derive the relation 

^-d\ = I T^dS. 
This relation and the above expression of Hi give 

Hi = j2 Res -^-dS. 

Q =i Pia A - Ti 

This leads to the second expression of the defining equation of the prepotential. 



32 



8.3 Proof of Proposition 

The idea of proof is the same as the standard one based on Riemann's bilinear relation. 
Integrability conditions to be checked are the following: 



dHj OH,, OHj <%/ dbj <% 7 



(8.4) 



dTi dTj' d ai &Tj d ai da/ 

Here bj, bj, Hi and Hj are used for simplifying notations; in the verification below, we 
have to return to the integral expressions of these quantities. 

These three relations can be verified in much the same way. We show the derivation 
of the second relation in detail. First, differentiating the integral formula 



by ai gives 



Hi = y — 1 -t= I -fi^dS 

3 V 27rv^T Jc ja A - Tj 



0Hj y.—^=f t%^- 

<-> . / T T„ \ TV 1 



da! a 27r v / ^T Jc ja A - Tj 
Similarly, from the integral formula of bi, 

By Riemann's bilinear relation, dflj and dujj satisfies the relation 

(*) == / VtjdujT = ( £ dVLj £ dujj — £ dVL,-. £ dour) . 

KJ 2n^lJdA 3 j\Ja k 3 Jb k Jb k 3 Ja k J 

Here A is a simply connected surface with boundary (4g-gon) obtained in a standard way 
by cutting the Riemann surface of the spectral curve along 2g paths. The boundary <9A is 
oriented in the direction encircling interior points anti-clockwise. We now evaluate both 
hand sides. As for the left hand side of (*), £lj has poles at the points Pji, • • • , Pj r in A 
whereas duji has of course no pole. Therefore the contour integral along dA splits into a 
sum of contour integrals along Cji, • • • , Cj r . Furthermore, the singular behavior of Qj in 
a neighborhood of P ja is such that 

Qj = —0j a / (A — Tj) + non-singular. 

Therefore 

Now consider the the right hand side of (*). The integrals of dVlj along A K ^s all vanish 
by the normalization conditions of dflj. Its integrals along B K 's do not vanish in general, 
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but now the integral of duj along A K 's vanish except for K — I; in the case of K = I, 
the latter integral is equal to 1 by the normalization conditions of duj. Thus 



RHS of (*) = -* dtlj. 

These relations imply the second integrability condition. 

In the verification of the first integrability condition, it is convenient to express dfli as 
a sum of meromorphic differentials dQ ia with just one pole at P,i a and normalized by the 
same condition of vanishing A-periods. The integrability condition can then be reduced 
to Riemann's bilinear relation for dVt ia and dVLjp. 

The third integrability condition is the easiest to verify. This case, however, is the 
most interesting in the context of the Seiberg-Witten theory. Differentiating bj now by 
aj gives 

— = f dui = T u . 

oaj JBj 

This is equal to dbi/daj, because the period matrix T is symmetric. (One can also 
directly deduce this conclusion by applying Riemann's bilinear relation to duj and duj. 
This is indeed a usual proof of the relation T LJ = Tjj !) 
This completes the proof of the proposition. 

As a corollary of the proposition and the final part of the proof, we obtain the following 
result: 

Corollary 2 The second a-derivatives of T — J-(T, a) coincide with the matrix elements 
of the period matrix: 

" =T U . (8.5) 



da j da j 

Thus our treatment of the prepotential incorporates all essential aspects of the prepoten- 
tial of Seiberg-Witten type. 

8.4 Final Remarks 

(i) The definition of prepotential shows a link with the notion of r function. Recall that 
Hi is equal to a logarithmic derivative of the r function of the Schlesinger equation. In 
the e-dependent formulation, 

Note, however, that log r and T are by no means in a simple proportional relation, because 
some averaging operation intervenes. 
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(ii) Unlike the prepotentials in topological conformal field theories and the Seiberg- 
Witten theory, our prepotential T seems to possess no manifest homogeneity of degree two. 
This issue seems to be linked with the unusual structure of dS that we pointed out in the 
end of the last section. Presumably our prepotential (and the modulation equation) will 
have hidden homogeneity, which emerges after introducing more deformation variables. 

9 Conclusion 

The body of this paper consists of two part. The first part is concerned with the derivation 
of the modulation equation. This equation is expected to describe slow dynamics of the 
spectral curve in isospectral approximation, in the sense of Gamier, of the (e-dependent) 
Schlesinger equation. Although our method for deriving this modulation equation is 
heuristic, we believe that this gives a correct answer (under a suitable choice of the 
symplectic homology basis). The second part of this paper is devoted to a complete 
description of solutions of the modulation equation, as well as the notion of prepotential. 
We have been able to obtain the following fundamental results: 

• General solutions are obtained by an inverse period map. The period map is given 
by period integrals of Seiberg-Witten type. 

• The modulation equation can coexist with a deformation equation of Seiberg-Witten 
type. They altogether form a commuting set of flows on the (^-dimensional moduli 
space of spectral curves. 

• A prepotential can be defined on this extended commuting flows. This construction 
is parallel to already known examples of prepotentials. 

These results, too, strongly support the validity of the modulation equation. 

A number of problems are left unanswered. The most crucial is of course the issue of 
validity of our derivation of modulation equation. We have presented a few fragmental 
ideas in this direction. Krichever's averaging method seems to be the most powerful 
and universal approach to this problem; we, however, have no idea how to calculate the 
averaged quantities to obtain period integrals. Presumably, Vereschagin's method will 
be a hint to this issue. Also, we would like to stress that generalizing the averaging 
calculation of Flaschka, Forest and McLaughlin to non-hyperelliptic spectral curves is 
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still a very important problem; the spectral Darboux coordinates of the Montreal group 
should be used in place of the classical auxiliary spectrum. 

The methods presented in this paper will be generalized to other types of isomon- 
odromic problems. An immediate, but also interesting generalization is the case with 
irregular singular points. For instance, if an irregular singular point of Poincare rank one 
is added to A = oo, the resulting isomonodromic problem becomes the so called JMMS 
This equations exhibits a remarkable "duality" |TT| , which will be inherited 



equation 



to the modulation equation. This issue will be discussed in a separate paper ||42|| . An 
even more interesting direction is a generalization to isomonodromic problems on an ellip- 
tic (and higher genus) Riemann surface. Examples of such isomonodromic problems have 
been constructed by Okamoto's group in a geometric framework [J43| . Recently, Levin and 
Olshanetsky presented another framework based on the method of Hamiltonian reduction, 



and pointed out a link with the Whitham equation 44 
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